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We study a quantum spherically symmetric object which is based on radial plasma oscillations. 
Such a plasmoid is supposed to exist in a dense plasma containing electrons, ions, and neutral 
particles. The method of creation and annihilation operators is applied to quantize the motion 
of charged particles in a self-consistent potential. We also study the effective interaction between 
oscillating particles owing to the exchange of a virtual acoustic wave, which is excited in the neutral 
component of plasma. It is shown that this interaction can be attractive and result in the formation 
of the ion condensate. We discuss possible applications of this phenomenon in astrophysical and 
terrestrial plasmas. 
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I. INTRODUCTION 



The studies of quantum effects in the dynamics of charged particles in dense plasmas is a rapidly developing branch 
of the modern plasma physics [l| . Numerous effects such as the behavior of quantum dots , the exciton dynamics Q 
etc, unusual for classical plasmas, were recently reported to exist when quantum dynamics is taken into account. Note 
that quantum effects in astrophysical plasmas are also important in the evolution of compact stars (see, e.g., Ref. 0]). 
However, as it was mentioned in Ref. still there is a lack of understanding how to correctly account for the 
quantum dynamics in plasma physics. In particular, a careful analysis of applicability of any of the approaches for 
the description of quantum plasmas should be made. 

In the present work we shall discuss the quantum dynamics of a certain class of plasma objects. We shall study 
a spherically symmetric plasmoid based on radial oscillations of charged particles. Note that such a configuration of 
plasma oscillations was previously analyzed in connection to the studies of the Langmuir waves collapse Q. We also 
mention that spherically symmetric plasma oscillations can be a theoretical model of stable natural plasmoids 0413 ■ 

We start in Sec. [TT] with the formulation of a model for a quantum spherical plasmoid based on radial plasma 
oscillations. Then, in Sec. IIII1 the motion of charged particles is quantized using the general ideas of Sec. |TTJ In 
Sec. HVI we carry out a quantization of acoustic waves which are excited in the neutral component of plasma existing 
in our system. We describe the effective interaction between charged particles owing to the exchange of an acoustic 
wave and analyze when this interaction can be attractive in Sec.fV] In Sec. lVII we study the possibility of the formation 
of bound states of oscillating charged particles in a spherical plasmoid. Finally, in Sec. IVIII some applications are 
considered. In Sec. I Villi we briefly summarize our results. 



II. MODEL OF A QUANTUM PLASMOID 



One of the main approaches for the quantitative studies of quantum effects in dense plasmas consists in accounting 
for the quantum pressure term in the Euler equation, i.e., for instance, the pressure becomes [111 ]. 

where m is the particle mass and n is the number density. Note that the modification (|2.1[) is likely to account for 
small corrections to the dynamics of classical plasmas. For example, there is an additional term in the dispersion 
relation for Langmuir plasma waves, uj 2 = w;?i(q) + fi 2 q 4 /4rrij (see, e.g., Ref. [l2|). where lj is the frequency of plasma 
oscillations, w c i(q) is the classical dispersion relation, m e is the electron mass, and q is the wave vector. Thus this 
approach for the description of the plasma dynamics seems to be the quasi-classical method for the studies of quantum 
effects in plasma physics. Therefore, to describe the dynamics of a strongly correlated system of charged particles, 
this method is unlikely to be valid and the new approach has to be developed. 

If we are interested in the studies of a quantum object based on oscillations of charged particles in plasma Q, the 
main problem in this description is the construction of a ground state of this many body system. The ground state 
based on the particles wave functions taken in the form of plane waves, as it is typically made in the condensed matter 
physics, is inappropriate since it does not fully reflect the dynamical features of the system. Thus we should build 
the ground state using the wave functions corresponding to the oscillatory motion of charged particles. 

The dynamics of plasma oscillations on the classical level was studied in Ref. [l3[ using the Lagrange variables 
method. In this formalism, the time evolution of a particular charged particle is analyzed. It was found that the 
separation of positive and negative charges in plasma results in the appearance of oscillations of lightest charged 



particles, i.e. electrons, with the Langmuir frequency cj e = y 4Tre 2 n^ /m e , where e > is the proton charge and n^ ' 
is the equilibrium number density of electrons. It means that the equation of motion of an electron has the form, 

f = -uj 2 e r. (2.2) 

Thus one can suggest that the perturbed motion of a single electron in plasma is governed by the quadratic potential 
V(r) — m e u> 2 r 2 /2. In this case one can reproduce the classical equation of motion (12.21) . 

We should also mention that the idealized picture described above when only electrons participate in plasma 
oscillations is never implemented. In a realistic case ions are also involved in the plasma motion. Although the 
frequency of ion oscillations, uji = u) e ^m e /mi, where m, is the ion mass, is significantly less than that of electrons, 
we should take into account the ion motion since quantum effects can be more strongly marked for ions. Note that 
accounting for the ion motion in the studies of plasma oscillations involves various electron-ion interactions. It makes 
Eq. (|2.2I) highly nonlinear. As a rule, these effects reduce the frequency of electron oscillations. Therefore we can 
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assume that electrons still perform oscillatory motion but with another frequency ujq < u e . Thus we will take that the 
motion of a charged particle in plasma is governed by the potential V(r) = rauj^r 2 /2, where m stays for the electron 
or ion mass and ujq is the oscillations frequency of charged particles. 

Let us assume that the dynamics of charged particles in a quantum plasmoid is also determined by the same 
potential V(r). Using this assumption, the construction of the ground state of this plasma structure will be made in 
Sec. IIIIl Of course, to have the correlated many particles system, the typical size of such an object should be quite 
small. We shall obtain some estimates of the length scale of these plasmoids in Sec IVIII Note that, in Ref. [Tj|, the 
cases of planar, cylindrical, and spherical plasma oscillations were studied. In our work we will be mainly concentrated 
on the description of spherically symmetric quantum plasmoids. 



III. QUANTIZATION OF THE CHARGED PARTICLES MOTION 



To describe the quantum dynamics of a particle in plasma we choose the Schrodinger picture, i.e. the operators 
of observables will be time independent. If we study the stationary states of particles, with the total wave function 
~ exp (— i^P) tpi where E is the particle energy, then the Schrodinger equation for the coordinate wave function tp has 
the form, 

E^ = Hf^, ^ fq = -|V + ^f^, (3.1) 

where we take into account that particles in a plasmoid move in the potential V^r), defined in Sec. ITT1 

Since the potential in Eq. (|3.1[) does not depend on angular coordinates, we shall look for the solution of this 
equation as ip = R(r)Yi m (9, <f>), where R(r) is the radial wave function and Yi m (6, <fr) are the spherical harmonics 
corresponding to the orbital and magnetic quantum numbers I and m. Introducing the new variable p = muj^r 2 jh 
and representing R in terms of the new function u, R(p) = e~ p / 2 p l / 2 u(p), we get the differential equation for the 
function it, 

d 2 u (Z , \ du ( E 3 l\ 
^ + {2 +l - p )dp- + {2^-l-2) U = - ^ 

The solution of Eq. (|3.2p has the form, 



- # (-5K + i + r' + 5-')' ( "» 

where $(a, b, z) is a confluent hypergeometric function. 

For the solution (|3.3I) to be finite we should require that the particle energy obeys the following relation: 

Era = fejo f 2n + i + |^ , (3.4) 

where n = 0,1,2,... is the radial quantum number. We will be mainly interested in the spherically symmetric 
solutions corresponding to I = 0. In this case the properly normalized total wave function has the form, 

, M 1 (-1)" /m 3 a; 3 \ 1/4 r mu a1 H 2n+1 (y^pr) 

where H n (z) is the Hermite polynomial, \a is the spin wave function, and a = ±1 is the spin variable. 
To secondly quantize the system, we introduce the operator valued wave functions, 

i>{r) = X! ®™ipna(r), (3.6) 
no- 
where a na - is the annihilation operator. One can also write down the analogous expression for the Hermiteian conjugate 
wave function which will contain the creation operator a' na . We shall suppose that ions are half integer spin particles, 
i.e. they are fermions like electrons. Thus, the operators a na and satisfy the canonical anticommutation relation 

= 5 nn '5 aa i for both electrons and ions. 
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Using Eqs. (13.11) and (I3.5|) we can express the total energy of the system of charged particles as 

E cp = J d 3 n/-t (r)H iq ^(r) = ^E^S^, (3.7) 



where 

E n = hu Un + , (3.8) 
is the particle energy which accounts for only the radial motion. Finally we replace the total energy 

E cp — > H cp — ^ E n a na a na , (3.9) 

no 

with the secondly quantized Hamiltonian. 

Now we can calculate the energy, Eq, of the ground state of the system. If the total number of particles, N, is even, 
using Eq. p.9p we get for Eq, 

E = 2 fwo ( 2" + g ) = 2?lw o ( n l + o 7lF + o ) ' ^ 3 - 10 ^ 
If AT is odd, we analogously obtain that 

E = 2 huj (in + + ftivo (^2n F + |) = 2ftw o ^| + ^n F + , (3.11) 

Note that, if N — > oo, both Eqs. (|3 . 10[) and p. lip approach the same limit Eq — > 2hoJon F . Here n-p stays for the 
maximal number of an occupied energy level. We can call it the Fermi number, in the analogy to the Fermi momentum. 

Let us define the effective size of a plasmoid, R, as the position of the last maximum of the function |^„ F (r)| 2 . If Tip is 
great, this maximum is approximately achieved at the classical turn point [141 ]. i.e. when Hloq (2tt,f + 3/2) = mtu^R 2 /2. 
Accounting for this fact, we get that 

mu R 2 

n F « (3.12) 

Using Eq. p,12j) we can express the ground state energy in terms of the plasmoid radius as Eq w m 2 u Q i R 4: /8h. 

Note that our plasmoid model is based on the linear approximation of noninteracting particles. As we mentioned 
in Sec. [TTl if we account for various (classical or quantum) electron-ion and electron-electron interactions, it will lead 
to the nonlinear terms in the equation of motion (|2.2I) of a particle and impose a certain length scale in the system. 
One the examples of such interaction is later discussed in Sec. [V] 

The ground state of a spherical plasmoid was constructed under the assumption of the negligible particle tempera- 
ture. The temperature corrections will wash out the lowest energy levels with E n < fceT, where fee is the Boltzmann 
constant and T is the electron or ion temperature. Taking into account the fact that particles possessing the lowest 
energy levels are concentrated near the center of a plasmoid, it means that plasma oscillations are the most intensive 
in a spherical shell. The internal radius of this shell is determined by the thermal effects and the external one by 
the nonlinear interactions between particles. Such kind of the spherical plasmoid structure is in agreement with the 
numerical simulations made in Ref. Ilal 



IV. QUANTIZATION OF THE ACOUSTIC FIELD 

As we mentioned in Sec. IIII1 the temperature in plasma should be not quite high for the existence of a spherical 
quantum plasma structure. It means that, along with the charged particles, there should be a quite significant 
fraction of neutral atoms or molecules in plasma. Plasma oscillations will inevitably result in collisions with neutral 
atoms leading to the generation of acoustic waves. In this Section we shall quantize acoustic waves inside a spherical 
plasmoid. 
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The hydrodynamics equations, which govern the acoustic waves propagation, have the form, 

^ + V(n„v„) = 0, 
at 

O -J 

-S5T + ( V » V K + V Pn = 0, (4.1) 

at n n m n 

where n n is the number density of neutral particles, v„ is the neutral gas velocity, p n {n n ) is the pressure of the gas, 
and m n is the neutral particle mass. Let us expand the parameters Ti rL1 v n , and p n as 

n n =n, ( i 0) + m H , 

v„ =vi H , 

Pn=Pi?+(^-) K-nW) + --, (4.2) 



dn r , 



where rin and are the equilibrium values of the number density and the pressure as well as n\ and vi arc the 
perturbations of the number density and the velocity. Accounting for Eq. (|4.2p , we can linearize Eq. (I4.1[) and present 
it in the form, 

^+ni°)(V-v 1 ) = 1 

-i + Jy Wi = °- ( 4 - 3 ) 

i In 

where 



m,. 



dp n 



1 / 2 V V dn n J n 

n V ^ / U 



(4.4) 



is the sound velocity. 

Now we can get the wave equations for vi and n\ from Eq. (|4.3|) . 



at 2 



c s 2 V(V • vi) - 0, (4.5) 



d n\ 



dt 2 

The spherically symmetric solution of Eq. f|4.6[) has the form, 



<V 2 n x = 0. (4.6) 



m(r,t) = e-^/fcW, AM = ^, (4.7) 

Ax 



where = c s fc. Note that the coordinate function obeys the following identities: 



r/*(r)/p(r) = * 8 (k-p), / ^fk(r)Mr')=d 3 (r-r'), (4.8) 



(27r) 3 -"" ' J (2tt) 

where delta-functions should be understood as <5 3 (k — k') = S(k — k')/Ank 2 . 

Following Ref. [la ], we introduce the potential of the velocity as vi = Vipi. Using Eqs. (|4.5[) - (|4.7[) we obtain the 
general expressions for ni and (pi as 

"'('•^/^(^)' ,2 ' i/ '(' , ^'"'' + ^ i "'')- 

^ = ^K^+^e^), (4.9) 



(2tt) 3 / 2 V2m„nW Wfc 
where 6^ and b\ are the Fourier coefficients. 
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The total energy of the acoustic field has the form, 

E ph = i [ d 3 r 



III,: 



(4.10) 



Now, if we replace bk and b\ in Eq. (|4.9[) by the annihilation and the creation operators, bk and b\, we can express 
i? p h, which should be regarded as the secondly quantized Hamiltonian of the phonon field, in the following way: 

E p h — > H p h = J d 3 k hujk b\bk + divergent terms, (4.11) 

To derive Eq. (|4. 1 1[) we suggest that the operators bk and b\ obey the canonical commutation relation for bosonic 
operators, [bk; b\i\- = J 3 (k — k'). Note the "divergent terms" in Eq. (|4.11[) . containing 5(0), can be removed by the 
normal ordering of operators. 

Note that, in this Section, we used the Heisenberg picture where operators turn out to be time dependent, cf. 

Eq. (14.91) . To obtain the operator in the Schrodinger picture, (bk) = &Sh, which will be used in Sec. [V] we should 

V / Sh 

make the canonical transformation, 

S Sh = oxp (-i^tj S H cxp (i^tj , (4.12) 

where &h = (bkj ■ Analogous transformation should be made for b\. Using Eq. (|4.11[) one can check that the 
transformation (|4.12p is equivalent to setting t = in Eq. (|4.9j) []~7j j . 



INTERACTION BETWEEN CHARGED PARTICLES AND ACOUSTIC WAVES 



We have already mentioned in Sec. IIVI that charged particles can interact with neutral atoms, which are present 
in the system. Thus charged particles should also interact with acoustic waves, which are generated by spherically 
symmetric plasma oscillations. In this Section we derive the secondly quantized Hamiltonian which describes this 
interaction. Then we exclude the acoustic part of the total Hamiltonian and reduce it to the nonlinear charged 
particles Hamiltonian. 

We suggest that charged particles scatter off the density perturbations of neutral atoms or molecules. Thus the 
energy of the interaction of an charged particle with an acoustic field has the form [l8| , 

V(r) = y dV ^(r-r')ni(r'), (5.1) 

where K(r — r') is the energy of the interaction between a charged particle placed at r and a neutral particle which is 
at r'. Using Eq. (I5.1[) we can derive the secondly quantized Hamiltonian of the interaction between charged particles 
and acoustic waves as 

^cp-ph = J d 3 rV3 f (r)U(r)^(r) = K J d 3 r^ (r)^(r)fti(r), (5.2) 

where ijj and h\ are given in Eqs. (|3.6[) and (j4.9[) respectively. To derive Eq. (|5.2I) we suggest that K(r — r') = 
KqS 3 (t — r'). This approximation corresponds to a contact interaction between charged and neutral particles. 
Using Eqs. (|3.6|) and (|4.9j) we can cast Hei-ph to the form, 

— — Dns ( fc )4er«s CT (bk + b\j , (5.3) 



where 



(0) 



1/2 



D ns (k) = K k ( ) / d 3 r ^„(r)i(f)/ t (r), (5.4) 
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is the matrix element. With help of Eqs. (13.51) and (|4.7[) we can compute D nm {k) at big n and s in the explicit form, 



K I nn^mujQ 
is) 1 / 4 y 2m n uj k ^ 

sign(2Vn-2v/i-C)-l], (5.5) 



1/2 



Dns(k) ^tt^iTI [sign(2^ - 2Vi + + sign(2v^ + 2^5 - £) 

4(ns)V 4 \ 2m n uJk I 



where £ = k^/h/mu! . 

To exclude the acoustic degrees of freedom we make the canonical transformation of the Hamiltonian in Eq. ([5 



Hep— ph ^ Hcp — ph 



e 



It can be seen that, after this transformation, the total Hamiltonian, H = Hq + /f cp _ p h, transforms into — > 
+ f [-ffcp-phj where -Ho = ^cp + -Hph- Averaging over the acoustic ground state and assuming that there are no 
external phonons, we get the total Hamiltonian in the form, 

7„ x 3 ^ D ns (k)D n > s >(k) 

7b (7 nn' ss' <j<j' 

Note that the Hamiltonian (|5.7p contains both the free charged particles term and the interaction between charged 
particles owing to the virtual acoustic wave exchange. Let us examine when this effective interaction can be attractive. 
We remind that charged particles in a plasmoid oscillate with a rather high frequency ~ ujq. In collisions with neutral 
atoms these charged particles will generate acoustic waves with the typical frequency ~ wo It means that the 

phonon energy in Eq. (|5.7[) is fiw^ ~ hujQ. Using Eq. (|3.8[) we that E n ~ fouj^n. It means that the effective interaction 
can be attractive only if E s = E„ or s = n, i.e. charged particles keep staying on the same energy levels. 



VI. BOUND STATES OF CHARGED PARTICLES 



In the previous Section we have found that the exchange of a virtual acoustic wave can result in the attractive 
interaction between oscillating charged particles in plasma. In this Section we study this process in details and show 
that under certain conditions charged particles can form a condensate. 

Before we proceed it should be noted that charged particles with parallel spins cannot approach one another close 
enough to effectively interact. It happens because of the Pauli principle. We remind that ions are also supposed to 
be fermions. Thus we should exclude the contribution of such charged particles from the effective interaction (15.71) . 
It mean that the summation in the nonlinear term should be made over the opposite spin indexes a' = —a. 

We should also mention that up to now we studied the case of constant number of particles. To avoid this restriction, 
we use the standard technique of shifting the energy levels E n — > e„ = E n — fi, where [i is the chemical potential of 
the system. Now the number of particles can change, but we should calculate the chemical potential. 

Finally, taking into account all these comments we can transform the Hamiltonian (|5 . T[) to the following form: 



- V F a f 



(6.1) 



where the new matrix element can be explicitly calculated on the basis of Eq. (|5.5p as 



Fnn '- Vo 7^ ' y °~^w!hrJ • (6 - 2) 

It is worth mentioning that F nn > > for any n and n'. It means that the effective interaction described by the 
Hamiltonian (|6.1|) is really attractive. 
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To diagonalize the Hamiltonian (16.11) we introduce the new operators Ao n and A\ n by means of the Bogolyubov 
transformation, 

a n+ = A 0n u„ + v n A\ n , a„_ = u n A ln - v n A\ n , (6.3) 

where the real coefficients Uk and Vk satisfy the relation + v\ = 1 . Note that the new operators obey the canonical 
anticommutation relations for fermion operators. Substituting these operators in Eq. (16.11) we get the Hamiltonian in 
the form, 

H = E' + H 2 + H' 2 + . . . , (6.4) 

where 

E' = 2j2<[e n -F nn ul], (6.5) 

n 

is the ground energy, 

i?2 = Y, [ e » K - V D + AF nnU 2 n vl] (A\ n A ln + ij„i „) , (6.6) 

n 

is the energy of the new quasiparticles, and 

H' 2 = 2 u n v n [e„ - F nn (u 2 n - vl)] (A ln A 0n + Al n A\^ , (6.7) 

n 

is the nondiagonal part of the quasiparticles interaction. Note that in Eq. (|6.4p we omit terms higher than quadratic 
ones. 

The nondiagonal interaction (|6.7p can be vanishing in two cases. Firstly, when either u n = 1 and v n — or 
u n = and v n = 1. It corresponds to the trivial solution, which is equivalent to the transition to the particles-holes 
representation. Note that the ground state energy (|6.5I) is not negative for this trivial solution: Eq > 0. 

The second situation is implemented if the attraction between charged particles is rather strong, i.e. when 

F nn >\e n \- (6.8) 
In this case we can choose the coefficients and Vk as 

°» = K 1+ i:)- ■i-K'-s;)- m 

The ground state energy (|6.5p now becomes negative, 

ffo = -E (F "9F 6n)2 <°- ( 6 - 10 ) 

It means that the state of the system, corresponding to the new quasiparticles, is more favorable energetically than 
in the trivial case discussed above. The diagonal part of the Hamiltonian (|6.6I) has the form, 



H 2 = F nn (A\ n A ln + 4„i n) • (6.H) 



Using Eqs. (|6.2p and (|6.1ip one can find the energy of quasiparticles E' n = 2VQ/y/n. 

The chemical potential of the system can be calculated using the expression N = (&' \N\Q ), where is the 
wave function of the ground state satisfying the conditions 

A „|^)=0, and i ln |^) = 0, (6.12) 

and 

N = J2^L&na, (6-13) 
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is the number of particles operator. 

Using Eqs. (|6.3p . (|6.9|) . (|6.12p . and (|6 . 1 3[) we can express the number of particles as 



^ = 2 E^= E 1 



(6.14) 



where the bound of the condensate in the energy space, rip, is defined by Eq. (|6.8j) . The quantity n' F is analogous to 
the Fermi number introduced in Sec. Mil Basing on Eqs. (|3.8|) and (|6.14j) we get the system of equations, 



N =— - H n ' 



15 



V ' 



2n F + -J - -JL, (6.15) 
which define the quantities n F and /it as functions of N. To derive Eq. (|6.15l) we use the approximate identity, 



E«° 

n<n' 



i a dn = 



(n F ) a+1 
a + 1 



(6.16) 



which is valid at big n F . 

The analysis of this Section shows that the exchange of a virtual acoustic wave between oscillating charged particles 
in a quantum spherical plasmoid results in the formation of their condensate. There is a pairing of particles with 
oppositely directed spins. These particles should be at the same energy level. Note that this new state of plasma 
has less energy of the ground state compared to the situation of unpaired particles. This process is analogous to the 
formation of Cooper pairs of electrons in a metal. We remind that, if the temperature of a metal is significantly low, 
two electrons can form a bound state owing to the exchange of a virtual phonon, which is a quanta of the vibration 
of a crystal lattice of a metal. This process underlies the phenomenon of superconductivity. 

Note that besides the oppositely directed spins, two electrons in a metal should have opposite momenta 20]. In 
our case the situation is analogous. Indeed, at big radial quantum numbers, n» 1, the asymptotic expansion of the 
charged particles wave function (13.5[) can be expressed as 



ipn(r) ~ - sin 
r 



muiQ-n \ 1 

r = 

h J 2ir 



mujon 
exp 2n I — : — r 



mu} n 
— exp — 2u I — : — r 



3.17) 



Eq. (|6.17l) means that any steady state of a charged particle in a spherical plasmoid is a superposition of converging 
and divergent spherical running waves. The pairing of two charged particles with opposite spins happens between the 
states corresponding to the running waves with opposite momenta. 

Let us demonstrate by means of the explicit calculation that the formation of a bound state of two charged particles 
is possible in our system. Suppose that there is the state of two charged particles having opposite spins, 



i*) = E c 



(6.18) 



where c„ 

a n <j\§o) 
is 



are the expansion coefficients and |$o) is the ground state corresponding to the filled lowest energy states: 
= 0, cf. Sec. Mil Using the Hamiltonian (I6.1[) we can show that the energy corresponding to the state (|6.18|) 



= (*|#|*) =2j2(En - F nn ) |c 



(6.19) 



where E n is given by Eq. (|3.8p . Thus the energy of the state (|6.18p is less than the sum of energies of two non- 
interacting charged particles, with F nn being the binding energy. 

Analyzing Eq. (|6.19p it is interesting to mention that, in contrast to the formation of Cooper pairs in metal, where 
bound states are formed near the Fermi sphere surface, in our case all the energy levels are subject to the particles 
condensation. This fact can be explained by the very convenient choice of the ground state wave functions (|3.5p . 
which maximally account for the dynamical features of the system (see also Sec. [IT])- Of course, we should remind that 
our analysis is valid for relatively big radial quantum numbers n. 
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VII. POSSIBLE APPLICATIONS 



To analyze the possibility of the charged particles condensation in a spherically symmetric plasma structure we 
can use Eq. (|6.19l) which gives one the energy of the bound state of two particles with opposite spins. Note that 
this expression is valid just before the condensation since it is based on the ground state |$q) rather than on |$q), 
defined by Eq. (|6.12[) . The condensation may happen if E n < F nn . Using Eqs. (|3.8[) and (16.21) . one obtains that this 
constraint is equivalent to 



Supposing that the Fermi number of a plasmoid before the condensation, given in Eq. (13. 12|) . is smaller than the 
constraint (17.11) . we get the following upper bound: 



on the radius of a plasma structure. 

Note the constraint (I7.2[) is very conservative. It means that charged particles on all energy levels can form bound 
states simultaneously. On the contrary, we may require that only particles at lower levels are subject to condensation 
and then this process spreads to higher energy states. Thus in a realistic case the upper bound of plasmoid radius (|7.2[) 
may be significantly changed towards its enhancement. However the analysis of the dynamics of this phase transition 
requires a separate special study. 

To evaluate the plasmoid radius one should define Kq, which does not depend on the macroscopic plasma charac- 
teristics. We remind that we approximated the potential K (r) by a delta-function. For this kind of potential, using 
the Born approximation, one can calculate the total cross section of the charged particles scattering off the neutral 
particles as a — m Kq/wH . Note that in our approximation the cross section does not depend on the particles energy, 
which is not the case for a realistic potential. 

Let us first examine the possibility of charged particles condensation inside a spherical plasmoid in a very dense 
medium corresponding to the inner crust of a neutron star (NS). The number density of neutral particles, i.e. neutrons, 

at the bottom of the inner crust is ~ 10 38 cm" 3 [2l|. We will be mainly interested in the studies of protons 
condensation. The number density of protons in the NS crust strongly depends on the equation of state of the NS 
matter. According to Ref. [22:] it can be about several per cent of the neutron density. In our estimates we take 

that rip ' 1 ~ 10 36 cm" 3 . Thus we get that at such a density the Langmuir frequency for protons is wo ~ 10 21 s _1 . 
The sound velocity also depends on the equation of state of the nuclear matter. Nevertheless we may take it as 
c s ~ 10 9 cm • s _1 [23|], which corresponds to the chosen density of the NS crust. 

Protons in the NS crust are highly degenerate. Their Fermi energy does not exceed several MeV. At such typical 
energies the cross section of the proton-neutron scattering is approximately constant and equals to 2 x 10~ 23 cm 2 24]. 
Thus we can also evaluate Kq. Finally, using Eq. (|7.2[) we get a plasmoid with R cr « 4.6 x 10" 12 cm is a subject of 
the protons condensation and possibly the superconductivity. The obtained value for the critical radius means that 
there can be several hundreds of energy levels inside t he p lasmoid. This result is in agreement with the fact that there 
can be a proton superconductivity in the NS matter [25j ]. 




Note that the energy corresponding to the zeroth level, 31twq/2, turns out to be of the order of the Fermi energy 
for protons. Thus thermal effects will not influence the plasmoid dynamics. 

It should be noticed that the cross section of the electron-neutron scattering is several orders of magnitude smaller 
than that of the proton- neutron scattering (26j . Moreover the Langmuir frequency for electrons is much higher than 
that for protons. It means that the critical volume of a plasmoid is too small to contain enough number of particles. 
Thus the electron superconductivity is unlikely to exist in a spherical plasmoid inside NS. 

Although one might expect that the charged particles condensation can happen only in plasmas with very high 
densities, which may be encountered only in astrophysical media, we may discuss the situation when this phenomenon 
occurs in terrestrial conditions. Let us discuss the case when a plasmoid appears in liquid water with ~ 10 23 cm -3 . 
We suggest that the number density of singly ionized water ions, is nf^ ~ 10 21 cm -3 . The Langmuir frequency for 
water ions is uio ~ 10 13 s _1 . The sound speed in water is c s ~ 10 5 cm- s _1 . Note that the total spin of a neutral water 
molecule is integer. Thus, if we study a singly ionized water ion, its spin should be half integer, i.e. it is a fermion. 

Unfortunately, the scattering of water ions on neutral molecules of water is not very well studied. Anyway we 
can assume that the total cross section of such a scattering cannot be less than that for electron scattering on 
water molecules which can be ~ 10" 14 cm 2 in the eV electron energy range [13]. Thus using Eq. (I7.2[) we get that 
R CI sa 4.6 x 10" 7 cm. Again we can see that there are several hundreds of energy levels inside the plasmoid. 




(7.1) 




(7.2) 




11 



As in case of a plasmoid existing in the crust of NS, the condensation of electrons in a spherical plasma structure 
involving water ions is unlikely to happen because of the very high frequency of electron oscillations. The fact that the 
exchange of a virtual acoustic wave results in the cohesion of ions rather than electrons was also noticed in Ref . [l8[ . 

The typical energy of a stochastic motion of an ion is fceli ~ 1CP 14 erg, where Ti m 300 K is the ion temperature. 
Using the above estimate for the oscillations frequency, we get that the energy of the zeroth energy level is also 
- 10~ 14 erg. Thus we can disdain the thermal effects in a spherical plasma structure with water ions. 

We have shown that the condensation of ions is possible in a plasmoid which is a radial oscillation of charged particles 
in plasma. It is also clear that ions play a subdominant role in plasma oscillations because of their low mobility. We 
may put forward a hypothesis that the condensation of electrons can also happen provided the frequency of their 
oscillations is significantly reduced owing to nonlinear electron-electron or electron-ion interactions. However, this 
issue requires an additional special study. 

The described condensation of water ions inside a plasmoid may result in the appearance of superconducting phase in 
plasma. Note that the idea that the plasma superconductivity may explain the stability of a natural plasma structure, 
called ball lightning (BL) (28[, was previously discussed in Refs. [29|,[30j]. As we have seen above, for the condensation 
of charged particles to happen, the typical diameter of a plasmoid should be very small, ~ (10~ 7 — 10 -6 ) cm. The 
plasma structures of the similar size, based on quantum oscillations of electrons, were studied in Refs. [§. [131 . It 
should be noted that previously the magnetic interaction (3l| and the exchange of a virtual Langmiur wave [32| were 
considered as possible mechanisms which underlie the plasma superconductivity. 

It should be noticed that the existence of the superconducting state of plasma inside BL can result in the presence 
of strong nondecaying electric currents in this kind of objects. It may explain the fact that sometimes the appearance 
of a natural BL magnetically affected heavy metallic objects like a church bell [33| . 

According to our estimates, the condensation of water ions can happen in a quite dense plasma with ~ 10 cm 
and nf^ oj 10 21 cm -3 , which corresponds to a liquid water. Thus our our results cannot be directly applied for the 
description of an atmospheric BL. We remind that the number density of saturated water vapor at the temperature 
~ 300 K is rin ~ 10 17 cm -3 , which is several orders of magnitude less than the value used in our work. Nevertheless 
our results may be used in the interpretation of the experiments with electric discharges in liquid water [U, [35[ , where 
luminous spherical objects, resembling a natural BL, were generated. 



VIII. CONCLUSION 



In conclusion we mention that in the present work we have developed the theory of quantum spherical plasmoids. 
Such a plasma structure is based on radially symmetric oscillations of charged particles in plasma. In contrast to the 
previous approaches to account for quantum effects in plasmas, which seem to be quasi-classical, our method involves 
quantum mechanical description of the electrons and ions evolution in a self-consistent potential which governs the 
oscillatory motion of charged particles. Thus our approach is valid for the description of dense quantum plasmas 
containing strongly correlated charged particles. 

In Sec. lIIII we have applied the general formalism to the quantization of charged particles in a spherically symmetric 
plasmoid. In particular we have derived the secondly quantized Hamiltonian (|3.9j) and have constructed the ground 
state of the system. In our work we have discussed the situation when, along with oscillating charged particles, 
there is a neutral component in plasma. In Sec. IIVI we have secondly quantized the field of acoustic waves which are 
inevitably excited inside a plasmoid. Then, in Sec. |Vj we have considered the interaction between charged particles 
and acoustic waves. We have derived the secondly quantized Hamiltonian of charged particles which accounts for 
the effective interaction owing to the exchange of virtual acoustic waves, cf. Eq. (|5.7|) . We have also considered the 
situation when this interaction can be attractive. In Sec. I VII we have discussed the attractive interaction between 
charged particles in details and have shown that it can result in the transformation of the ground state of the system. 
We have demonstrated that charged particles, which are supposed to be 1/2-spin fermions, tend to form singlet bound 
states which are more favorable energetically. This phenomenon is analogous to the formation of Cooper pairs in a 
metal which results in its superconductivity. 

In Sec. IVIII we have discussed the possible applications of the condensation of charged particles in plasma owing to 
the exchange of virtual acoustic waves. Firstly, we have considered the case of a very dense plasma corresponding to 
the inner crust of NS. It has been shown that the proton component of plasma can be subject of the condensation. 
This result is in agreement with the previous findings that proton superconductivity may well happen in the NS 
matter. Secondly, we have studied the possibility of condensation of singly ionized water ions inside a plasmoid in 
the terrestrial conditions. We have shown that this phenomenon can happen provided a plasma structure is created 
in a liquid water and the density of ions is quite high. We have also considered the implication of our results to the 
theoretical description of stable natural plasma objects as well as to the resent experiments where plasmoids were 
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generated in electric discharged in water. 

We have obtained that the condensation is unlikely to occur in the electron component of plasma. It happens 
because of the quite high frequency of electron oscillations. Thus electrons will excite acoustic waves ineffectively. 
Moreover we should take into account the higher electron temperature compared to that of ions. It means that the 
formation of bound states of electrons will be washed out by their thermal motion. This result is consistent with 
the claim of the authors of Ref. 18] who found that the exchange of an acoustic wave leads to the cohesion of ions 
rather than electrons. We also notice that our mechanism of the attractive interaction between ions is different from 
that recently proposed in Ref. (3(| , where the shielding of the repulsive potential of an ion by the electron cloud was 
considered. 
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